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Resume 

This article aims to find explicit congruences between Dirichlet characters and 
gives various results on how to find some effectively on a computer, ft ends with 
concrete examples putting those ideas in application. 

1 Introduction 

Le but de ce texte est, a partir du choix de deux entiers N ^ 2 et M ^ 2, d'obtenir 
une congruence modulo m pour les caracteres de Dirichlet modulo N. 

Plus precisement, si G est le groupe des caracteres de Dirichlet modulo TV, alors on 
cherche des entiers algebriques explicites (a x ) x6 c tels que : 

S %, V] a x x{x) = mod M 

XGG 

La motivation derriere ce travail est le theoreme suivant, que Ton trouve dans la 
these [5] de l'auteur : 

Theoreme 1 On suppose que I' on dispose d'une combinaison lineaire finie de ca- 
racteres de conducteurs Mg, qui verifie une congruence : 



^7xX = OmodMf0 

x 

alors : 

CI E Txff (/. a ' S > 3, x)L Ui3 , (j + l,x)=0 mod M^ h ^O 



x 

oil : 



H(f,a,S,j,x) 
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1 



Ce theoreme abstrait sur lequel on ne souhaite pas s'etendre (il y a beaucoup de 
notations et de concepts qui sont loin du sujet du present article), affirme l'existence 
de certaines congruences pour des valeurs speciales de fonctions arithmetiques (des 
fonctions L de formes modulaires, tordues par des caracteres), qui sont verifiees si Ton 
dispose de congruences pour les caracteres. On souhaite done construire des exemples 
concrets de ces dernieres, pour un jour verifier experimentalement dans quelle mesure 
les congruences entre les valeurs speciales sont optimales. 

Le plan de ce travail est le suivant : 

- on commence par expliquer comment ramener le probleme initial a un calcul de 
noyau matriciel ; 

- on explique ensuite comment determiner le noyau, avec le logiciel SAGE0), d'une 
fagon completement automatique d'une part (mais tres inefficace) et de facon plus 
guidee d'autre part (mais tres manucllc) ; 

on determine des congruences dans des cas particuliers explicites ; 

- finalement, on donne le code source des fonctions utilisees. 

Je souhaite remercier Alexei Pantchichkine pour son soutien indefectible, et William 
Stein pour m'avoir donne acces aux machines du reseau math. Washington. edu0. 

2 Reduction a un calcul de noyau matriciel 

On se donne comme dans l'introduction un premier entier N ^ 2, et on considere le 
groupe G des caracteres de Dirichlet modulo N, que Ton liste sous la forme xi, ■ ■ ■ ,Xn 
oh n = \G\. Si on definit m = N — 1, on peut alors representer les valeurs prises par 
tous les caracteres simultanement sous la forme d'une matrice de taille (m, n), avec : 

= Xj(i - 1) 

De cette fagon, si (ai, . . . , a n ) est une famille de scalaires et x £ [0, N — lj, la quan- 
tity Y^i=i a iXi( x ) es t l e coefficient en ligne x + 1 dans le vecteur colonne A(a± . . . a n ) T . 

Notons que cette matrice initiale iaun noyau trivial : e'est le resultat bien connu 
d'independance des caracteres. 

Cependant, comme les caracteres considcrcs sont a valeurs dans les racincs de 
l'unite, si on fixe F un corps de nombre qui les contient, la matrice A est a coeffi- 
cients dans l'anneau des entiers O de ce corps. Pour des raisons pratiques, on choisira 
bien sur le corps le plus petit possible, done un corps cyclotomique de plus petit degre 
possible (que Ton notera d). 

Si on se donne maintenant un entier M ^ 2, la recherche de congruences modulo 
cet entier a coefficients dans O revient a chercher un vecteur-colonne V T = (ai . . . a n ) 
tel que AV £ MO n . Si on note B l'image de A dans l'anneau des matrices a coefficients 
dans 0/(M), on est ramene a chercher les elements du noyau de B 7 puis a les relever. 

1. version 5.4, disponible sur http//www . sagemath. org 

2. Les calculs presentes dans cet article sont tous aises et rapides sur une machine simple, mais les 
tatonnements et experiences necessaires a leur recherche ont parfois necessite l'acces a une puissance de 
calcul plus consequente ; les machines en question ont ete financees par "National Science Foundation 
Grant No. DMS-0821725" . 



Le probleme des congruences pour les caracteres de Dirichlet est ainsi ramene a 
un probleme lineaire de calcul de noyau de matrice a coefficients dans un quotient 
d'anneau d'entiers. 

3 Calcul automatique du noyau 

D'apres un resultat classique (voir par exemple le lemme 6 de Particle de Birch [1]), 
l'anneau O est le S-module engendre par les puissances d'une racine primitive de 
l'unite. L'anneau Oj (M) est done un module de type fini sur S/MS, done un ensemble 
fini : il suffit de considerer tous les vecteurs possibles pour obtenir le noyau. 

C'est theoriquement tres satisfaisant, mais en pratique : TLjMTL est de cardinal M, 
done 0/(M) est de cardinal M d , et les vecteurs a considerer sont done M dn . 

Pour le cas le plus simple que Ton presentera de fagon detaillee en 15. 2\ N = 5 ct 
M = 16, cela represente 4294967296 vecteurs a considerer. Comme de plus les calculs 
vectoriels dans l'anneau en question sont en plus relativement lents, il est clair que 
cette approche est d'une utilite assez limitee. 

Donnons malgre tout le code source effectuant le calcul ; il utilise une fonction 
matrix_of_Dirichlet_group discutee plus loin : 

def integer_mod_iter(d,M,zeta) : 

for coeffs in product (*tee (range (M) ,d) ) : 
result = 
for ii in range(d): 

result=result+coef f s [ii] *zeta**ii 
yield result 

N = 5 
M = 16 

A = matrix_of _Dirichlet_group (N) 

= A.base_ring() . ring_of _integers () 

d = A.base_ring() .degreeO 

Q = D . quotient_ring(0 . ideal (M) , 'a') 

zeta = Q.gens()[l] 

B = matrix(Q,A) 

m,n = B . dimensions () 

print (M**(d*n)) 

for X in product(*tee(integer_mod_iter(d,M,zeta) , n)): 
img = B*vector(Q,X) 
if img. is_zero() : 
print (X) 



4 Calcul guide du noyau 



4.1 Principe general 

Le probleme du calcul du noyau d'une matrice a coefficients dans un corps est 
simple ; on peut par exemple penser a un passage sous forme de matrice echelonnee 
reduite (c'est, par exemple car on le trouve partout, l'algoritlime 7.3 dans le livre de 
Stein [3])- On peut aussi songer a la notion de diviseur elementaire, mais elle necessite 
un anneau principal et le calcul effectif un anneau euclidien. 

La matrice qui nous interesse est a coefficients dans un anneau quotient, qui n'a 
pas de tres bonnes proprietes : les fonctions habituelles et les algorithmes implementes 
dans SAGE ne fonctionnent done pas. 

L'idee va etre d'obtenir une relation de la forme B * R = L * E, ou L et R sont 
des matrices inversibles, et E est aussi proche que possible d'une matrice echelonnee 
reduite. Le calcul de vecteurs dans le noyau de E fournit alors des vecteurs dans le 
noyau de B, tout simplement en considerant leur image par R. 

Plus precisement, on va tenter d'obtenir E sous forme de blocs : 



ou le bloc est un bloc identite (done carre), et Q est un bloc (a priori rectangu- 
laire) dans lequel on garantit : 

- aucun coefficient n'est inversible ; 

- aucune ligne n'est nulle ; 

- aucune colonne n'est nulle. 

C'est le mieux que l'on puisse demander a SAGE de faire automatiquement ; c'est 
la raison pour laquelle la suite de la recherche necessite une intervention manuelle 
(discutee plus loin). 

Le calcul du noyau est alors ramene au calcul du noyau de cette matrice reduite E ; 
il fait intervenir : 

- des elements triviaux, qui correspondent aux colonnes nulles a droite (s'il y en 



- des elements non triviaux, qui correspondent au noyau du bloc Q (que l'on voudra 
done reduire a la main apres avoir obtenu une premiere version) . 

On appelera l'entier r (taille de la matrice identite dans la decomposition) le pseudo- 
rang : il correspond a un nombre de colonnes que l'on sait etre independantes sur 
l'anneau considere. II depend evidemment de la decomposition obtenue, mais s'il est 
egal a n, on est assure qu'il n'existe pas de congruences modulo M. Le nombre de 
colonnes nulles sera appele le noyau garanti. 

4.2 Algorithme de reduction automatique 

L'idee est d'amenager l'algorithme du pivot de Gauss, en l'analysant en termes 
d'operations elementaires sur les lignes et les colonnes. 




a) ; 



Si est un couple d'indices distincts, on peut definir Pjj la matrice de permu- 
tation associe a la transposition de i et j, et si a est un scalaire, on peut considerer 
Tij(a) la matrice de transvection qui est une matrice carree identite modifiee avec le 
coefficient a cn ligne i et colonne j, et si a 0, la matrice de dilatation Di(a), matrice 
identite modifiee avec le coefficient a en ligne i et colonne i. Les operations elementaires 
sur les lignes se decrivent par multiplication a gauche par ces matrices, et les operations 
elementaires sur les colonnes par multiplication a droite. Attention, ces matrices sont 
carrees, mais on ne fait pas apparaitre leur taille dans la notation. 

Le principe va etre le suivant : si on a reussi a ecrire B*R = L*E avec E quelconque, 
alors on a aussi les egalites suivantes, pour les operations sur les lignes : 



qui montrent que si une operation elementaire permet de « simplifier E >, on peut le 
faire quitte a modifier L (pour les operations sur les lignes, et par une operation sur 
les colonnes) ou R (pour les operations sur les colonnes, et par la meme operation sur 
les colonnes). 

Expliquons comment on travaillc : 

1. On part tout d'abord avec le triplet (L,E,R) = (I m ,B,I n ), qui verifie par 
construction B * I n — I m * E ; e'est l'invariant B*R = L*Ede l'algorithme. On 
commence avec k — 1 ; 

2. Si on trouve un coefficient inversible dans E en position (I, c), alors : 

(a) on permute les lignes I et k de E, et les colonnes I et k de L (ce qui conserve 
l'invariant) ; 

(b) on permute les colonnes c et k de E de R (pour l'invariance) ; 

(c) le coefficient unite se retrouve en (k, k) ; par dilatation, on divise la ligne k 
de E et on multiplie la colonne k de L par ce coefficient (meme remarque) ; 

(d) le coefficient est maintenant egal a 1 ; on ajoute dans E, avec multiplication 
par un coefficient la k-eme ligne puis la k-eme colonne aux autres lignes et 
colonnes pour que la ligne k et la colonne k ne contiennent que ce 1 (e'est 
un pivot) ; operations sur E que Ton compense sur L et R comme discute 
ci-dessus, toujours pour conserver l'invariant ; 

(e) on incrcmcntc k et on recommence (dans la mesure ou k $5 min{m,n}). 

3. Si on ne trouve plus de coefficient inversible, on se contente de faire migrer par 
des permutations les lignes et les colonnes nulles de B vers le bas et la droite 
(toujours en compensant sur L et R). 



B*R = LPij * PijE 

= LT itj (-a)*T itj (a)E 
= LDi(l/a)*Di(a)E 



et sur les colonnes : 



B * RP ttJ 
B*RT hJ (a) 
B * RDi(a) 



L * EPij 
L * ETij (a) 
L * EDi(a) 



La premiere partie de pivot de l'algorithme garantit que Ton a une allure 
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sans element inversible hors du premier bloc, et la migration finale des lignes et colonncs 
nulles garantit que le bloc central n'a plus ni ligne ni colonne nulle et fait apparaitre 
les parties droites et basses nulles. 

4.3 Discussion de ['implementation 

4.3.1 matrix_of_Dirichlet_group 

Cette fonction recoit l'entier TV ^ 2 choisi, puis calcule et renvoie la matrice associee 
au groupe des caracteres de Dirichlet, qu'il faudra ensuite reduire modulo un entier 
M ^ 2. C'est la seule qui ne soit pas contenue dans l'implementation de la classe. 

4.3.2 dirty _cached_is_unit 

La version du logiciel SAGE utilisee ne dispose pas d'un moyen de test de l'inver- 
sibilite dans un anneau quotient de l'anneau des entiers d'un corps de nombres ; mais 
propose une methode d'inversion sous reserve d'existence (avec levee d'exception si ce 
n'est pas le cas). On a done implements une fonction de test a l'aide de l'inversion. 

C'est tres couteux, mais aura le merite d'etre aise a remplacer lors devolutions 
futures du logiciel. Pour gagner un peu de temps, on memorise les resultats. Experi- 
mentalement, on gagne relativement peu avec cette mise en cache (quelques pourcents). 

Ellc retourne un booleen. 

4.3.3 Classe invariant_triplet 

On travaille comme on l'a vu avec des triplets de matrices (L,E,R). La presence 
de L n'a aucun interet a priori vu l'objectif de calcul d'un noyau, mais on la conserve 
tout de meme, pour deux raisons : 

- la premiere est qu'elle coute peu a maintenir done ne constitue pas une reelle 
gene, 

- et la seconde parce que le maintien de l'invariant donne un moyen de verification 
intercssant. En particulier, durant le developpement, de nombreuses coquilles ont 
ete reperees par des tests de l'invariant sur des matrices tirees au hasard. 

Ensuite, comme on le verra plus loin, l'algorithme automatique ne donne pas de 
tres bons resultats : une fois choisi le couple (M,N), il faut guider la machine. C'est 
la raison pour laquelle on a modelise la situation a l'aide d'une classe invariant_triplct, 
dont on initialise les instances avec la matrice B, et qui va se charger de gerer le triplet. 

Detaillons l'interface programmatique de cette classe : 

Initialisation La fonction d'initialisation regoit une matrice de taille quelconque, et 
initialise l'objet avec le triplet (L, E, R) — (I m ,B, I n ). 

Variables On peut inspecter les variables L, R et E. En particulier, l'inspection de 
E permet de decider quelles operations de reduction on va demander. 



Methode assert_invariant Verifie que l'objet respecte toujours l'invariant - extremement 
utile durant le developpement. 

Methode check_kernel_vector Regoit une liste de coefficients, qu'elle utilise pour 
creer un vecteur-colonne. Si ce vecteur colonne est dans le noyau de E, alors la 
methode retourne une paire constitute du booleen True et de l'image du vecteur- 
colonne par la matrice R (done le vecteur explicite donnant la congruence 
cherchee !) ; sinon elle retourne une paire constitute du booleen False et de l'image 
du vecteur-colonne par la matrice E (ce qui permet de voir a quel point on s'est 
trompe pour corriger). 

Methodes do_row_addition et do_column_addition regoivent deux indices i et j 
et un coefficient a. Elles realisent Taction de Tij(a) sur le triplet. 

Methodes do_swap_rows et do_swap .columns recoivent deux indices i et j, et 
realisent Taction de Ptj sur le triplet. 

Methode do_pivot realise la partie pivot de Talgorithme de reduction presente en !4.21 
page|U 

Methode do_migrate_zeros realise la partie de deplacement des lignes et des co- 

lonnes nulles vers le bas et la droite dans Talgorithme presente en 14.21 
Methode do_normalize appelle successivement les fonctions do_pivot et do_migrate jeros. 

5 Resultats obtenus 

5.1 Statistiques 

La reduction automatique, passant en revue les cas oil M £ [2, 20] et N e [2, 20]] 
(soit 361 paires), a trouve 64 cas oil le noyau est garanti comme reduit a zero (matrice I r 
occupant toute la largeur de la matrice), et oil done aucune congruence n'est possible. 

Les frequences des rangs garantis parmi les 297 cas restants est listee dans la tabled] 
et les frequences des noyaux garantis dans la table [2j On constate done que les matrices 
rencontrees meme avec de petits parametres ne se reduisent pas bien automatiquement : 
la recherche de congruences necessite beaucoup d'afhnage manuel. 

Table 1 - Frequences des pseudo-rangs 



pseudo-rang 


frequence 


1 


180 


2 


54 


3 


27 


6 


27 


7 


9 



5.2 Congruences pour N = 5 et M = 16 

On choisit de se concentrer sur le cas N = 5 et M = 16 car e'est un cas de 
congruence de degre 4 le long de 2, et la dimension est petite : cela semble done un bon 



Table 2 - Frequences de noyaux 



noyau garanti 


frequence 





279 


1 


13 


3 


5 



choix dc premier exemple. Remarquons que cet exemple minimaliste pour la methodc 
guidee est hors de portee de la recherche exhaustive, comme on l'a vu en [3] 
Voici le debut de la session de calcul : 



N = 5 
M = 16 

A = matrix_of _Dirichlet_group(N) 

= A.base_ring() . ring_of _integers () 

m,n = A.dimensionsO 

Q = . quotient_ring(0 . ideal (M) , 'a') 

B = invariant_triplet(matrix(Q, A)) 

B . do_normalize () 

Que Ton poursuit ensuite a la main par : 

zeta4=Q.gens() [1] 
B . do_row_addition(3 ,2, 1) 
B . do_row_addition(3 , 1 , -1) 
B . do_column_addition(3 ,1,1) 
B . do_row_addit ion (1,2,-2) 
B . do_row_addition(l ,3,1) 

A partir de la, la matrice E est assez simplifiee pour que l'on devine des vecteurs 
du noyau ; apres avoir elimine les doublons, il reste les vecteurs : 

B . check_kernel_vector ([0,0,0,8]) 
B . check_kernel_vector ([0,0,4,0]) 
B . check_kernel_vector ([0,8,0, 4*zeta4-4] ) 

On trouve done des vecteurs de congruences explicites : 
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5.3 Congruences pour N = 7 et M = 15 

Voici le debut de la session de calcul : 



N = 7 
M = 15 

A = matrix_of _Dirichlet_group(N) 

= A.base_ring() . ring_of _integers () 

m,n = A. dimensions () 

Q = D.quotient_ring(0.ideal(M) , 'a') 

B = invariant_triplet(matrix(Q, A)) 

B . do_normalize () 

(a parte : M dn vaut ici 129746337890625) 

Que Ton poursuit ensuite a la main ; au bout de quelques etapes, on est ramene a 



E 



et 



On trouve done des vecteurs de congruences explicites : 
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6 Code 

def matrix_of _Dirichlet_group (M) : 
G=Dir ichletGroup (M) 
characters = G.listO 
field = G.base_ring() 

return matrix (field, M, len(characters) , lambda lig.col: characters [col] (lig) ) 



cached_units = O 

def dirty_cached_is_unit (a) : 

if cached_units .has_key(a) : 



return cached_units [a] 
else : 

result = False 
= a.base_ring() 
try: 

b = 0(l/a) 
result = True 
except : 
pass 

cached_units [a] =result 
return result 

class invariant_triplet : 

def __init__(self , B) : 
self.B = B 
self .E = copy(B) 
self .m, self .n = B . dimensions () 
self.L = identity_matrix(B .base_ring() , self.m) 
self.R = identity_matrix(B .base_ring() , self.n) 
self . assert_invariant () 

def assert_invariant (self ) : 

if self .B*self .R != self .L*self .E: 
raise ArithmeticError 

def check_kernel_vector(self , coef f s) : 

v = matrix(self .E.base_ring() , self.n, 1, coeffs) 
Ev = self .E*v 
if Ev. is_zero() : 

return (True, self.R*v) 
else : 

return (False, Ev) 

def do_row_addition(self , i , j , a) : 

self . E . add_mult iple_of .row (i , j , a) 
self . L . add_multiple_of _column( j , i , -a) 

def do_column_addition(self , i , j ,a) : 

self . E . add_multiple_of _column(i , j , a) 
self .R. add_multiple_of _column(i , j , a) 

def do_swap_rows (self , i , j ) : 
self . E . swap_rows (i , j ) 
self . L . swap_columns ( i , j ) 



def do_swap_columns (self , i , j ) : 
self . E . swap_columns ( i , j ) 
self . R . swap_columns ( i , j ) 

def do_pivot (self ) : 
k = 

while k < min(self .m,self .n) : 
found_pivot = False 
1 = k 
c = k 

while 1 < self .m and c < self .n and not found_pivot: 
if dirty_cached_is_unit (self . E [1 , c] ) : 

found_pivot = True 
if not found_pivot: 

1 = 1+1 

if 1 == self.m: 
1 = k 
c = c+1 
if found_pivot: 

self . do_swap_rows (1 , k) 
self . do_swap_columns (c ,k) 
coef f =self .E[k,k] 
self .E.rescale_row(k, 1/coef f ) 
self .L.rescale_col (k, coef f) 
for 1 in range (self .m) : 

if 1 != k and self .E[l,k] != 0: 
coeff=self .E[l,k] 
self . do_row_addit ion (1 ,k, -coef f) 
for c in range (self . n) : 

if c != k and self .E[k,c] != 0: 
coeff = self.E[k,c] 
self . do_column_addit ion (c,k, -coef f) 

k = k+1 
else : 

break 

def do_migrate_zeros (self ) : 
target = self.m-1 
considered = target 
while considered >= 0: 

if self . E [considered, :] .is_zero() : 
if considered < target: 

self .E. swap_rows (considered, target) 
self . L . swap_columns (considered, target) 
target = target-1 
considered = considered-1 



target = self.n-1 
considered = target 
while considered >= 0: 

if self . E [ : , considered] . is_zero () : 
if considered < target: 

self .E. swap_columns (considered, target) 
self . R. swap_columns (considered, target) 
target = target-1 
considered = considered-1 

def do_normalize (self ) : 
self . do_pivot () 
self . do_migrate_zeros () 
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